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Abstract. We propose an experimental setup to measure the work performed in a 
normal-metal/insulator/superconducting (NIS) junction, subjected to a voltage change 
and in contact with a thermal bath. We compute the performed work and argue that 
the associated heat release can be measured experimentally. Our results are based on 
an equivalence between the dynamics of the NIS junction and that of an assembly of 
two-level systems subjected to a circularly polarised field, for which we can determine 
the work-characteristic function exactly. The average work dissipated by the NIS 
junction, as well as its fluctuations, are determined. From the work characteristic 
function, we also compute the work probability-distribution and show that it does not 
have a Gaussian character. Our results allow for a direct experimental test of the 
Crooks-Tasaki fluctuation relation. 
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Fluctuation relations are central to our present understanding of statistical mechanics. 
Their long and distinguished history goes back to at least the studies by Callen and 
Welton [T], Green |2] and Kubo |3], which were inspired by the works of Einstein on 
the Brownian movement |1] and of Johnson [5] and Nyquist [6] on noise in electrical 
circuits. 

The derivation by Jarzynski of a rather general non-equilibrium work relation [7], 
linking the average of the exponential of the work being performed on a system with 
the equilibrium free energy difference between initial and hnal equilibrium states of the 
system, is of particular interest. Subsequently, Crooks mi obtained the Jarzynski 
equality from a relation between the probability of a given amount of work being 
performed on a system and the probability that the system performs the same amount 
of work on its surroundings if the time-reversed protocol is undertaken. Such heightened 
interest in non-equilibrium work relations has not only been fuelled by novel theoretical 
advances in out-of-equilibrium dynamics, but also by experimental progress in preparing 
and probing non-equilibrium evolution (see liniEIGSIGSlIIl] for reviews on the subject). 

In the present paper, we propose a relatively simple yet realistic experiment 
based on a proposal by Crooks [15], which would allow for a direct test of the 
Crooks-Tasaki fluctuation relation. The experimental system we consider consists of 
a normal-metal/insulator/superconducting (NIS) junction between a superconductor 
and a normal metal, which is initially short-circuited and is subjected to a given 
voltage protocol (see below). We establish that for the proposed protocol, the full work 
distribution-function has non-Gaussian character with a non-standard decay exponent. 
We also compute the hrst two moments of such distribution, which can be determined 
experimentally by measuring the average heat released and its variance. 

In order to obtain a closed expression for the characteristic function of the work 
distribution, we derive an equivalence between the dynamics of an NIS junction and 
the one of an assembly of two-level systems subjected to a circularly polarized held. 
Using this mapping, the work distribution is determined. For realistic parameters of 
experimentally available NIS junctions as e.g. those of Ref. [151, our hndings show 
that at cryogenic temperatures, even the tiny amounts of heat released will result in a 
measurable volume change for a probe connected to the junction. 

The structure of the paper is as follows: Sect. places the recent developments 
briehy outlined above into their proper historical context and reviews the theoretical 
tools necessary to study work-huctuation relations at the quantum level. In Sect. we 
introduce the Hamiltonian of the NIS junction and describe in detail the work protocol 
applied to the system. In Sect. we use results that are derived in Appendix A and 


Appendix B to compute both the average work dissipated as well as its huctuations and 


the full work distribution. [Appendix C contains an extension of the results presented 
in [Appendix Bj and Appendix D contains a discussion of a method of measurement of 
the tunnelling amplitude that describes the NIS junction. Section contains the main 
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result of our paper, where we compute the numerical value of the average heat released 
on a NIS junction such as the ones discussed in Ref. [16]. We leave to Appendix E the 


discussion of the experimental techniques that can be used to measure such a quantity. 
In particular, we show that the heat that is released can be measured by determining 
the volume change of a probe that absorbs it. We provide estimates of the relative 
volume change, which should be measurable using capacitive methods HZIIIH]. Finally, 
in Sect. we present our conclusions. For increased readability, most technical details 
have been relegated to the appendices. 


2. Perspective and previous works 

In this section, we give a brief historical overview of the evolution of the subject of non¬ 
equilibrium fluctuation theorems. Besides the early contributions already mentioned 
above, a statement of a general fluctuation theorem involving the free energy of a gas 
of hard-spheres, interacting at large distances via the Lennard-Jones potential, can be 
found in the work of Zwanzig [19]. A subsequent important development was the work 
of Bochkov and Kuzovlev [20l [21], which also concerned the relation of probabilities 
between a given work protocol and its time reversed version, but which adopted a 
different perspective to that of Jarzynski and Crooks, namely the former authors 
considered the generalised force that performs work on the system as being external 
to the dynamical description of the system, rather than internal. Another signihcant 
contribution was the work of Evans and Searles (see [22] and references therein), which 
considers the relation between the probability for a dynamical trajectory characterised 
by a certain value of a dissipation function and its time-reversed counterpart. Regarding 
developments concerning the validity and extension of the Jarzynski-Crooks relations at 
the classical level, see the experimental works [23l [211 EH EHl [271 EHl EHl [30l EU [32l [33] 
and the theoretical ones [SH EH EH EH EHl |39] , the above list not being exhaustive. 

A quantum generalisation of the Crooks equality was obtained by Tasaki m and 
also by Kurchan [H] (for the case of cyclic protocols). A version of the Jarzynski 
equality valid for quantum systems was also derived by Yukawa [12], albeit treating 
work at a quantum level as an observable, a view that was shown not to be correct 
[IH] - Such a view lead in the past to a debate concerning the validity of such fluctuation 
relations at the quantum level [IH EH EG]- The extension of the Jarzynski equality 
to isolated systems, using the dehnition of the work distribution function of Tasaki 
and Kurchan, is due to Mukamel mi- A quantum generalisation of the work of 
Bochkov and Kuzovlev referred above, as well as a clarihcation of the relation of 
their work to that of Jarzynski and Crooks, can be found in a paper by Campisi et 
al. [IH]- Extensions of non-equilibrium work relations to isolated systems (micro- 
canonical fluctuation theorems) are given in [IH [50], whereas work relations valid 
for arbitrary open quantum systems were obtained in [5l]. The last reference shows 
the validity of such theorems for systems that can arbitrarily exchange heat with 
the bath while a work protocol is being applied to them (see also below). Other 
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Figure 1. (a) NIS junction, (b) Schematic representation of the voltage protocol 
considered 

theoretical developments concerning quantum systems can be found in the references 

[SaEaliaiEslESlETllMlESllmlEIlEllEallHlEsllsniETlEHlEglEQllZIlEi, where 
again the list is not exhaustive. 

At the experimental level, several proposals to demonstrate the validity of the 
Crooks-Tasaki relation in the quantum domain have been presented, including the 
use of a series of projective measurements 1731 , the measurement of this relation in 
the optical spectra of systems subjected to weak quenches m, or the use of qubit 
interferometry [Z5l ESI EZ!. The former as well as the latter proposal have been 
successfully implemented, see [ZSl E9]. An example of the experimental conhrmation 
of a generalised version of the Jarzynski equality in the quantum domain was discussed 
in EH. Finally, a discussion of possible solid state experiments performed on quantum 
heat engines is given in [72] . 

The above mentioned proposal by Crooks [I5] relies on an indirect method of 
observation namely, the measurement of the heat released by a quantum system after 
a cyclic work protocol has been performed on the system. A specihc set-up to measure 
the heat released by a circuit including a resistor and a Cooper-pair box, based on the 
temperature increase of a local probe, was already presented in [HO], where the Cooper- 
pair box acts as a single two-level system. The present proposal has the additional 
advantage that the heat released scales with the contact area between the normal-metal 
and the superconducting dims. 

The work probability distribution, P /3 (IT, r), is dehned, for a class of work protocols 
where the system is decoupled from the thermal bath during the application of the 
protocol (see below), of duration r, as jlO] 

Pf(W,r) = ^ E I (n.T |C/(t,0)| m, 0 ) ^ , (1) 

where | m, 0) are the eigenstates of the system’s Hamiltonian at the beginning of the 
protocol, P(0) (which is given for the NIS junction by (|^ at t = 0, see below), and 

are the corresponding eigenvalues, and where | n, r) are the eigenstates of the 
system’s Hamiltonian at the end of the protocol, H{t) (which is given by the same 
equation at t = r), with Ef being the corresponding eigenvalues, with ft = l/^ksT) and 
Z{ft) = Tr(e“^'^^°^). The operator P(r, 0) is the time-evolution operator of the system 
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during the protocol, i.e. in the interval [0,r], which can be written as a time-ordered 
exponential of the full Hamiltonian of the system, given for a NIS junction by eq. ([^. 

The work characteristic-function of the NIS junction is the Fourier transform of eq. 
Q, and is given by 




Tr 


W (r, 0) f/(r, 0) p(/?) , 


( 2 ) 


where p{(3) = Z{f3) is the density matrix of the junction before the application 

of the protocol. 

The unitary evolution of the system in the interval [0, r] reflects the work protocol’s 
restriction mentioned above, i.e. the NIS junction is disconnected from the thermal bath 
at f = 0 and remains adiabatically insulated while the protocol is being applied, being 
reconnected to the thermal bath at t = r and undergoing equilibration through the 
exchange of heat with the bath afterwards. It is possible to lift such a restriction on the 
protocol by explicitly considering the interaction of the system with the thermal bath 
[SI], but we will not consider such an extension here, see however the discussion below. 

The knowledge of the work characteristic-function allows one to compute all work 
moments. If one performs the analytic continuation v = ifd in the above definition and 
uses the cyclic character of the trace operation, one obtains 

T) = ZAI3)IZ(I}) = , (3) 

where Zr{fd) = Tr{e~^^^'^'^) is the partition function for a junction in equilibrium, 
described by the Hamiltonian H (r), AF being the difference in free-energies between the 
final equilibrium state (after thermalisation at a time t > r) and the initial equilibrium 
state at t = 0. Since in our case the Hamiltonian goes through a cycle, AF = 0. 
Expressing through its Fourier transform, we obtain 

(g-/3W) _ _ g-/3AF _ ^ ^ ^4^ 

J —oo 

which is a statement of the Jarzynski equality, which the work probability distribution 
satisfies by definition. Considering the relation between Q^iv^r) and Gfsiv + i(d,T), 
using again the cyclic character of the trace, it is possible to derive the identity 
F^(-IH,-r) = e-^^P^(hF,r) where by P^(hF, —r) we denote the work probability 
distribution for the time-reversed protocol. This is a statement of the Crooks-Tasaki 
relation. 

Note that in the case of cyclic protocols, one also has AU = 0, where U denotes 
the internal energy of the system. Therefore, from the first law of thermodynamics, 
AU = {W)p — AQ, where AQ denotes the average heat released by the system on 
which the work protocol is applied, we obtain AQ = (hF)^, which shows that for cyclic 
protocols the measurement of the released heat allows for the determination of (W)^, 
as stated. 
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The NIS junction is composed of a normal-metal film, a thin insulating layer and 
a superconducting him. We consider that such a junction can be described by an 
Hamiltonian consisting of three terms: a Fermi-gas describing the normal-metal him, a 
BCS superconductor describing the superconducting him and a tunnelling Hamiltonian 
coupling the two [HI [82l [83l [HI [85l |86l [87] , which gives rise to induced superconductivity 
in the normal-metal him through the proximity ehect (see [HE] and references therein). 
We use the Hamiltonian 


H{t) = H^{t) + Hs + X. (5) 

The operator Hn{t) = (^k 0(t)) cl^ka^nka is a free-fermion Hamiltonian describing 

the normal-metal, with 0(t) being the voltage diherence across the junction, whose value 
changes in time. The operator 

Hs = ^ cJka-Cskcr + ^ ('^kCsktCs_k4, + ^kCs-kJ^C^kt) ) (6) 

k,(T k 

is the Hamiltonian describing the BCS superconductor, where the pairing function Ak 
is diherent from zero on a band of width 2huD around the Fermi-level, where un is the 
superconducting material’s Debye frequency. Finally, the third term of Eq. ([^ is given 
by 

X = ^ r{k) (c^kaCnka + 4-k-aCs-k-^) , (7) 

k,(T 


and represents the tunnelling process across the junction. We consider the overall 
tunnelling matrix element to be small with respect to the gap function Ak. Such matrix 
element is taken to be invariant under time-reversal, i.e. = t~^{—k) [H]. The 

fermions’ momentum is directed along the plane of the films and we assume it to be 
conserved in a tunnelling process. We will measure the kinetic energy of electrons in the 
normal-metal him and in the superconducting him with respect to their common Fermi- 
energy /i. The super- or subscripts n and s refer to operators or quantities pertaining 
to the normal-metal or to the superconductor, respectively. The magnitude of the 
tunnelling amplitude will depend on the width of the aluminium oxide layer and its 
value can be adequately measured for a given device using the method discussed in 


Appendix D In this method, the power dissipated by the junction is measured at a 


constant applied voltage equal in value to the energy gap of the superconductor (divided 
by e), which is much larger than the value used in our protocol. 

We consider the NIS to be initially {t < 0) in equilibrium with an heat bath at 
temperature T and that the voltage 0 across the junction is zero. Such voltage is 
changed from 0 to a hnite value fmax and returned to zero within a hnite time-interval 
starting at t = 0 and ending at t = r. During such an interval the junction is decoupled 
from the heat bath. For simplicity, we consider that the voltage protocol is symmetric 
with respect to the mid-point of the interval t = t 12. However, the calculation presented 
in Appendix A is valid for more general protocols. Also, (fit) varies smoothly within a 
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time-interval t' and is returned to zero at the end of the protocol within the same 
time frame (see £gure[^(b)). 

Physically, for small values of the tunnelling amplitude, the term given by eq. 
([^ leads to the opening of a gap between the energy bands of the normal metal, of 
magnitude A"" A, leaving the energy bands of the superconductor largely unaffected. 
This is the proximity effect [88]. In our protocol, we place ourselves within a particular 
adiabatic limit, i.e. we will consider that the switching-time r' S> h/A, where A is 
the value of the gap-function within the superconductor. Thus, the voltage difference 
across the junction is held approximately constant within the microscopic time-scale 
of the superconductor and the change of voltage does not induce transitions which 
involve the superconducting bands. However, we also take r 3> h/A” 3> r' S> h/A. 
This implies that the dynamics is fully diabatic with respect to the separation between 
the two energy-bands in the normal metal and thus that the voltage protocol induces 
transitions between these bands, leading to work being performed on the system, see 
Appendix B 

Using the Hamiltonian for an NIS junction and the voltage protocol specihed, the 
work probability distribution can be obtained using Eq. ([^, as will be shown in the 
next section. 


4. Work dissipated due to the application of the voltage protocol to the 
junction 


In this section, we establish that the dynamics of the NIS junction is equivalent to 
that of an assembly of quantum spins in a time-dependent magnetic held. This entails 
an equality between the work probability distributions of both systems, from which it 
follows that one can compute the average heat released by the NIS junction, as well as 
its higher order huctuations. 

Starting from Eq. ([^, the explicit dependence of the Hamiltonian on the time- 
dependent bias voltage can be eliminated by a gauge transformation at the expense 
of acquiring a time-dependent hopping term between the normal and superconducting 
Elms. Thus, the Hamiltonian becomes 

H,{t) = HM + Hs + X{t) ( 8 ) 


with 

X{t) = t’ik) , (9) 

k,(j 

where $(t) = /q du(f){u). Anticipating that, in general, t^{k) <C A, one can, 
in the specihc adiabatic limit discussed above, perform a generalised Schrieffer- 
Wolff transformation that takes into account the explicit time-dependence of X{t). 
Furthermore, such a transformation can only be performed if eUnax ^ ^5 see [Appendix 


[A| After such a transformation, the metal and superconducting Elms efiectively decouple 
and only the transformed metallic excitations are subjected to the time-dependent bias. 
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d 


k,cr 




cl_u I -l- A. 


k4. 




( 10 ) 


where 


AJ = - 


|i'’(k)l"(«E + 5J 


((Sr-(Ei)^ ' 

iti(k)p + iti(k)p: 
21 (5E)= - (Eiy 1 


( 11 ) 


( 12 ) 


are, respectively, the renormalized kinetic energy of electrons and the induced 
superconducting parameter, in the normal metal, due to the proximity effect, = 
\J (^k)^ + being the energy of the superconducting excitations. This result is 

shown in detail in Appendix A[ We thus conclude that for spin-independent tunnelling 
matrix elements, the Hamiltonian describing the normal-metal is equivalent, in each 
subspace, to the Hamiltonian of a two-level system under the action 
of a circularly polarised held. Furthermore, the Hamiltonian that describes the 


superconductor does not contribute to the work-characteristic function, see [Appendix 

One can obtain the work characteristic function for a NIS junction simply by 
considering the product of individual characteristic-functions for two-level systems 


subjected to a circularly polarised held, see Appendix B The logarithm of the overall 
work characteristic-function, which is the generating function of the connected work- 
moments, is given by 

W/3(u,r) = ln^^(u,r) = { (1 “Pk) 


cosh 


Pk ■ 


i/3 + 2iv) VlAl^P + (f-)2 


cosh /3V|AjJ|2 + (^")2 


|A"lHeHnax)" 


with pk = 

calculation can be found in [Appendix B 


sm 


ASP + (5 


eVr 


max 


(13) 


The detailed 


The average work dissipated per atom in the normal-metal him, w, is obtained from 


the derivative of Eq. (13), evaluated at u = 0, 

(2 


w = 


2(eK 


max; 


AT. 


at 


tanh (/^Vl^kP + (^k)^ 

k 


X 


lAkl"' sin' I { v|AE|s + (a - eVmax)' 


|AE|s + (®'(|Ai;p + (5 


(14) 


where Afat is the total number of atoms in the normal-metal him. From Eq. 0 , one 
has that for weak-coupling between the two hlms, .^k ~ ■Ck- Moreover, if, as in the 
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system experimentally studied in [16], the dispersion relation is the same in both the 
normal-metal and the superconducting him (since the normal-metal him is composed of 
the same material as that of the superconducting one, but weakly doped with another 
metal), we have from Eq. (12) that . Assuming that both the order parameter 

Ak and |t(k)p only depend on k through their dependence on and are only non-zero 
in a vicinity of width Hud around the Fermi level of the material composing the hlms, 
where is fhe Debye frequency of the said material, one can transform the above 
summation into an integral and write the above equation as 

w = 2(el4iax)^ [ dfp{f) tanh f/5Y^|A”(0p + ^' 

J-hiJo V ^ 

(s vAhbkAAAAp) 


|A”(OP 


sm 


X 


yiMOF+e (|A”K)P + K - eV„,ax)^) 


( 16 ) 


where we have introduced p{f) = ^ ~ ^k)) fhe density of states per atom (and 

per spin-species) in the normal metal him. 


Diherentiating Eq. (13) twice with respect to v, at v = 0, one obtains the mean- 


square deviation and thus the relative deviation y/dw^/w of the work performed by the 
NIS junction during the protocol. If one transforms such result into an integral using the 


same assumptions we have used to obtain Eq. (15), this integral can also be computed 


numerically. In a similar fashion, higher moments can be calculated. 

The full characterisation of the work huctuations requires the determination of the 
work probability distribution, through the computation of the inverse Fourier transform 
of Writing the probability distribution in terms of the logarithm of the 


characteristic function, as dehned in Eq. (|13|), one has 
P^iW,T) = 


r+°^ dv ^-iWv+WB[v,T) 

—oo 27r 


( 16 ) 


Since W is an extensive quantity, i.e. it is proportional to A/'at, the integral can 
be computed using a saddle-point approximation. Minimising the argument of the 


exponential in Eq. (16), one obtains the following saddle-point equation 


f'hujD 


w = 


f — TlIjJd 


dip{C)d^S{f,x) , 


where w = W/Afat is the work per atom of the sample and S{f, x) is given by 


S{i,x) = ln<^ (1 -p(0) +P(0 ■ 


cosh 


(P) + 2x)^|A”(OP + e 


cosh 


/3j|A-(op+e 


(17) 


(18) 


where p{f) = ^ 


|A’^(OP(eVmax)" 


sm 


(AA”K)l" + (fNote 


that the solutions, x{w), of the saddle-point equation ( |l7| correspond to the analytic 
continuation of the exponent of the integral given in Eq. ( p!6| to imaginary values 
of n = —ix{w). Numerically solving Eq. (0 for a set of values of w (with the 
same assumptions as in the two calculations previously performed) and substituting 
its solution in Eq. (16), one obtains the result plotted in hgure for the logarithm 
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of PjiiW^T), divided by Nat- This calculation indicates that the work probability 
distribution is a stretched exponential for both positive and negative W. In the inset, 
the exponent associated with such a stretched exponential is determined for positive W 
(blue triangles). 

We have also subtracted the factor fiw from the logarithm of the probability 
distribution at negative W (green squares in the inset) and the function so obtained 
coincides with the logarithm of the probability distribution for positive IT, thus showing 
that the distribution obeys the Crooks-Tasaki relation at the level of the numerics. 



Figure 2. Logarithm of the work per atom distribution for a AlMn/AbOs/Al junction 

at r = 1 a:. 


The above hgure is a clear indication of the non-Gaussian nature of the work 
distribution, which may allow the reconstruction of the distribution above from 
experimental data, using the method that we propose in [Appendix E 


5. Calculation of the dissipated work in the NIS junction 

Our primary motivation for the present study was the analysis of the validity of the 
Crooks-Tasaki relation in an experimentally relevant setting. As discussed above and in 
Appendix E[ this can be accomplished by measuring the heat released in the NIS junction 


during a full cycle of the protocol. Since such measurements require a high degree of 
control of the energy flows occurring in an experimental set-up, it is advantageous that 
the released heat scales with the contact area between the superconducting and normal 
metal dims. 

We now consider the computation of the dissipated work released as heat by the 
NIS junction. As a concrete example, we will consider the junction studied in reference 
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[I6] . where manganese doped aluminium (the Mn concentration being 0.4% in mass) 
acts as the normal-metal electrode in the junction and pure aluminium acts as the 
superconducting electrode (with aluminium-oxide being the insulating material). 

In aluminium, the density of states is approximately constant in the vicinity of the 
Fermi surface yielding p ^ O.OllTcab^ per atom (and per spin) [HHl EO], in units of the 
Debye frequency, which is approximately ud = 10^^ Hz in aluminium. This metal is 
a BCS superconductor with a transition temperature, T^. = 1.20K [HI], and the gap- 
function given by A = 1.764 fc^Tc, where one uses the universal (BCS) relation [92], 
valid in the weak-coupling limit. For an intermediate coupling superconductor the gap is 
larger than predicted by the BCS relation. This is of no concern here, since we assumed 
r 3> r' 3> h/A ~ 3.6 x 10“^^ s. 

For concreteness, we consider the example of a protocol of duration r = 10“®s 
with an applied voltage Imax = 8 x 10“® V at a temperature of T = 1 K. Performing 
the integral in (15), assuming a typical tunnelling element t{k) ~ 0.05A, we obtain an 
average work w ~ 3.2 x lO'^^ha;^) per atom which amounts to approximately 0.7nJ 
of dissipated work per gram of aluminium. For these parameters, we also obtain a 
standard deviation of 5w‘^ ~ 8.4 x 10“^^(ha;£))^ corresponding to a relative deviation, 
VduP / w ~ 2.9 X 10^, which is a large value for such a quantity. 

The measurement of the released heat can be performed using several techniques, 
of which the measurement of the volume change of the him in isothermal conditions 
seems to be the most promising one, see Appendix E for a detailed discussion. 

A word is due with respect to the approximation that the dynamics is unitary during 
the application of the protocol, i.e. that no heat is exchanged with the thermal bath 
in the time interval t G [0,r], occurring only at later stages. The exchange of energy 
between the electronic degrees of freedom and the heat bath occurs mainly through 
electron-phonon inelastic collisions. The dependence of resistivity on temperature due 
to electron-phonon collisions is given, in metals, by p(T) = po(T/ToY, where Tq is 
a temperature of the order of the Debye temperature, which is 428 K in aluminium 
[93] . Using the value po = 3.875 x 10“®Dm, with Tq = 400K [M] for aluminium, and 
assuming a Drude like dependence of the resistivity on the electron-phonon relaxation 
time Tpe, p{T) = where m* = 1.10is the electron effective mass in 

aluminium l95 




^^ and n = 2.1 x 10^^ m”'^ is the electronic density in aluminium, we 
obtain rpe(T = 1 K) 5 x 10“^s, which is much larger than the protocol duration r. 
Hence, the said approximation is fully justified. Moreover, such relaxation time is of the 
order of the equilibration time that is necessary to wait after the end of the protocol 
before the volume deviation discussed in [Appendix E| can be measured. 


6. Conclusions 

Fluctuation-dissipation relations play a central role in thermodynamics and the 
approach to equilibrium, but the direct experimental verification of the various relations 
can often be difficult. In this work, we discussed a simple but realistic experimental set- 
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up that will allow to test the Crooks-Tasaki fluctuation relation. Our proposal consists 
of an NIS junction of a superconductor and a normal metal where the bias voltage is 
altered according to a specified protocol, see Fig. [^(b). 

As shown, the released heat scales with the contact area between the 
superconducting and normal metal Aims. Our proposal is thus amenable to experimental 
verification, e.g. by coupling the junction to a probe that isothermally absorbs the 
released heat and changes its volume as a result. Order of magnitude estimates of this 


effect for the aluminium-based junctions of Ref. [I6] were presented in Appendix E 


The calculation of the full work distribution function and its first moments clearly 
establishes an exponential behavior of such quantity that displays a non-Gaussian 
character with a non-standard decay exponent. 

Our results are based on an equivalence between the dynamics of a NIS junction 
and the dynamics of an assembly of two-level systems subjected to a circularly polarized 
held. This equivalence holds in a particular adiabatic limit in which the switching-on/off 
of the potential difference across the junction takes place in time scales much larger than 
the microscopic time scale of the system, determined by the value of the superconducting 
gap. 
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Appendix A. Mapping of the work characteristic-function to that of an 
assembly of two-level systems 

Following [96] and starting from the Hamiltonian ([^, we change representation, so as 
to write the original time-evolution operator f/(f, 0), determined by such Hamiltonian, 
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in terms of the time-evolution operator in a new representation, V{t, 0), with 0 < t < r, 

U{t, 0) = V{t, 0), (A.l) 

where <h(t) = f^ducjiu) and Nn = Sk.o- ^nka^niicr is the number of particles on the 
normal side of the junction. This change of representation corresponds to a gauge 
transformation that eliminates the dependence on 0(f) of Hn{t) and transfers such 
dependence to the tunnelling part of the Hamiltonian. The transformed operator, 
Hi{t) = jg given by Hi{t) = Hq + A A’(f), where A is a small 

dimensionless parameter that sets the scale of the tunnelling matrix element and which 
is explicitly written here for convenience of calculation. The operator Hq = Hn{0) + Hg 
describes the normal-metal him and the superconducting him at f = 0, in the absence 
of the tunnelling operator (7). On the other hand, the operator T’(f) is given by eq. ([^. 

The operator V{t, 0) can be expressed in terms of a time-ordered exponential 
involving the transformed Hamiltonian Hi{t). This exponential can be written in terms 
of an ordered product of exponentials, evaluated at increasing times, as 

N 


H(f,0) = nexp 


i=i 


it 


hN 




(A.2) 


where N ^ oo and ui = I ■ t/N. We now introduce a second-unitary transformation, 
dependent on both A and on t, by dehning a transformed Hamiltonian operator as 

% it) = % it) , (A.3) 

where the Hermitian operator S{t) will be adequately chosen below. Substituting this 


result in (A.2), we obtain 
t>(f,0) = e 


N r 

- exp 

1=1 ^ 


it 


hN 


H2{ui) 




JXS(O) 


. (A.4) 


Using the Baker-Campbell-Hausdorf formula, one can write, to order A^, and to order 
1/N 


^iXS{ui) -i\S{ui-i) ^ iX{S{ui)-S{ui-i))+^ [S{ui),S{ui-i)]- 


iX 

N 


dS I iX 
du^ 2 




(A.5) 


where [ , ]_ is the commutator of two operators and where the derivative of S{u) is 
computed at ui. Substituting (A.5) in (A.4), we obtain to order A^ and order 1/N, the 
result 


N 


v{t,0) = Yl exp 


i=i 


it 


hN 


n{ui) 


„iXS(0) 


(A.6) 


where, to order A^, the operator H(t) is given by 

?{(t) ^ (^0 + A N(t)) _xh^-'-^ 

(Lb ^ 
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2 

2 


s{t),-n^ + x{t) + i[s{t),Ho 


S{t ), X{t) 


(A.7) 


We choose the hermitian operator S{f) such that the linear term in A vanishes in 

(A.8) 


(A.7), i.e. such operator obeys the first order ordinary differential equation 


dS 


h— = X{t)+i[S{t),H^ 
whereas the transformed Hamiltonian reduces to 
fL{t) = Ho + ^-^[S{t),X{t) 


(A.9) 


One still needs to specify an initial condition for (A.8) to be properly defined, which 


we will do below. The above choice of S{t) corresponds to a time-dependent Schrieffer- 
Wolff transformation |9H], which is reminiscent of the work of 1100] . One can 
thus write the operator f7(t, 0) as f7(t, 0) = g-*AS(t) g) where 


r (t, 0) = T exp 


— - / duTLiu) 
h Jo 


(A. 10) 


which is valid to order and where 'H(m) is given by (A.9). Substituting this expression 


for U{t, 0) in (|^ and using the cyclic character of the trace operation, one obtains 
Qfi{v,T) = Tr 


0 ) ^-i\s(t) 


■ r (r, 0) p{(j) 


(A.ll) 


Under the action of, respectively, 5(0) and 5(r), H{0) and Hi{t) do not transform 
into 'H(O) and Tiir), as the factors that involve the derivative of S{t) at these points. 


and which appear in (A.7), are not present in the transformation. However, we may 


still fix the derivative of S{t) at one of these points, by the choice of an appropriate 

= 0. The explicit computation of Qjsiv^r) now 


initial condition. We thus choose ^ 


t=o 


requires the solution of (A.8), subjected to the condition that we have imposed at t = 0. 


In order to solve this equation, we write S{t) in the form 
Sit) = cLaCnka + Kit) ^kaC.ka ] 


kcr 


5Z['®k(^) ^s-k-a^n-k-a + ^kit) Cn-k-aCs-k-a \ i 

kcr 


(A. 12) 


where r^{t) and s^{t) are complex-valued functions, to be determined. Substituting 


(A. 12) in (A.8) and equating the terms pertaining to the same pairs of operators, we 


obtain the coupled system of first-order ordinary differential equations 

jjk — ie^{t)lh i (pn _ ps\ a , ^ A it 

df-h^ + ft Ski ’^k + ft ^k Sk 


dsl 

dt 


= |Ki; + ®sS+S<^AkrJ 


(A.13) 
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dr^ 

with the initial conditions 


(A.13) at t = 0, yield for r^(0^ 


dt 


t=0 

and 5^(0^ 




= 0. These conditions, once substituted in 
t=o 

, the result 


(A. 14) 


4(0) = 


it'^(k) (7 Ak 


This choice of parameters corresponds to the time-independent Schrieffer-Wolff 
transformation introduced in mu and such transformation can be applied to the 
Hamiltonian of the system at f < 0, since 0(f) = 0 in this interval. 

Writing the pairing function in terms of an amplitude and a phase, i.e. Ak = 
lAkle**^, where (f is independent of k, we define the new variables _ = 


u^{t) e aV’, y^{t) = v^{t)e Substituting these expressions in eq. (A.13) 


and differentiating the result, we obtain two decoupled second-order equations for x^{t) 
and y^(t), with appropriate boundary conditions. These equations can be easily solved. 

The solutions of such equations, expressed in terms of the original functions r^{t), 
Sk(f), are given by 


r^(f) = iT^(k) 

1 A 

^ hE? Jo 






(Eir 


cos{ult) 


Et 


sin(ci;kf) 


du ( e(j){u) - ^k - ^k ) 


and 


= 


. ?k“) sin(c(;k(f — u)) 


ia AkE(k) e^^k* 
1 


(A.15) 


- (^k 


sW 


cos(a;kf) 


Et 


sin(a;kf) 


hEf^ 


sin(a;^(f - «)) 


(A.16) 


Now, one can show, using the know properties of the exponential function and 
integration by parts, the following identity 


1 


hEt Jo 


du sin(a;^(f-M)) = 


- (^k)^ 




- |^cos(n;^f) - ^ • sin(n;^f)^ 

- ^ / dti e0(M) 

n Jo 

■ cos(a;k(f — u)) (A.17) 

- ^ ■ sin(n;^(f - u))^ . 

Using (A.17) in (A.15) and (A.16), one can write the solution of the system given in 

r(k) 


(A.13) as 


= ^E{k) (^k+ ^k) . fm I 

^ - {Eir 


- (E^y] Jo 


du e(j){u) 
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— I 


.. ?k(?k + «k) + |Ai 


Ei 


and 


s-m - - - 

- TTwvi 'iTPs\2 ^ 


■ sin(n;^(t - u))j 
r(k) cr Ak 


(A.18) 




{Kf 


mi? 


X “)) . j cos(a;k(t — u)) 


iK?] 


due(f){u) (A. 19) 


Et 


sin(a;k(t — u)) 


We stated above that we are considering protocols in which the applied voltage difference 
(fit) changes from 0 to Wnax in a time interval r' i.e. (fit) = Vuiax?(f/r'), where fj^x) 
changes from 0 to 1 in an interval of length 1. The integrals which enter the second 


term of both (A.18) and (A.19) are all, after the change of variables u —)■ u/rf of the 
form 


cldnax'r^ \ ^Ahl3i2£ZLvi/(u') -iiLtcr 

- - / duifiu)e fi '- >e ^ 

h Jo 


±E^)u 


(A.20) 


where T(m) = Jff dvif{v). The term if{u)e is a bonnded fnnction, whereas 

the exponential is at least eqnal to which in the adiabatic regime 

t' 3> h/A oscillates rapidly. Thns, provided that t ^ t' and eldnax'^V^ A 1, then 


(A.20) is zero by the Riemann-Lebesgne lemma. From these conditions, we also obtain 


elmax ^ A, i.e. the tnnnelling voltage applied shonld be well within the gap. Therefore, 


we conclnde that in the adiabatic limit and for t r', (A.18) and (A.19) rednce to 


r£(t) 


sm 




(A.21) 


it'^(k) cr Ak 


e n 


fm 


From snch a discnssion and since (f{T) = 0, it also follows that ^ 


and hence we obtain ^ 

at 


= 0. In snch a case, we can write eq. (A.11) as 


t=T 




dsZ 

r'^ k 

dt 


t=T 


Gi3{v,t) = Tr 


^(r, 0) 0) e-m)v 


(A.22) 


where 'H(r) is given by (A.9) (with t snbstitnted by r in that eqnation). Snbstitnting 
the operat or S {t) by its expression as given in (A.12), with r^it) and Sk(t) given by 

lere the operators TLg and 7in{t), 


(A.21 


, in (A.9), we obtain that 'H(t) = 'Hs + 'H„(t), w' 


pertaining, respectively, to the snpercondnctor and the normal-metal, are given by 


tts — ^ 'fko- ^Ika^ska + ^ A Ajj.Cs_l4Cskt) 5 


k,(T 


with 


= ^k - 




s\2 


Ak = Ak- 1 


- (^k 

W |tt(k)|2 + |tl(k)f 

2 


i^l? 


(El? 


(A.23) 

(A.24) 

(A.25) 































The distribution of work performed on a NIS junction 17 

being the renormalised kinetic energy and pairing function in the superconductor, and 

«»(() = i: C 


k,(T 


with 




+ E (a£ fifcty+a;; 

k 

A2|r(k)p(ej^ + a) 




.n_ Ak(|tt(k)p + |t^(k)p: 

“ - y - 


(A.26) 


(A.27) 


(A.28) 


(e£)2 - 

being the renormalised kinetic energy and induced pairing function in the normal-metal, 
due to the proximity effect. However, such an induced pairing function is multiplied by 
a time-dependent phase, due to the applied voltage across the junction. If one assumes 
that the tunnelling square amplitudes are independent of spin, i.e. |f^(k)p = |t'*'(k)p, 
the modihed kinetic energies are also independent of the spin one can, introducing in 
each ((k, j'), (—k, j,)) subspace the operators, 

K = cLfCnkt + “ 1 , (A.29) 

^k = ^Ikt^n-k^. 5 (A.30) 

kJ.Cnki' ) (A.31) 

whose algebra is isomorphic to the spin 1/2 algebra, write Unit), up to a constant factor, 
as 


^n(t) = E (^'k K + K d+ + A” dk ) 


(A.32) 


We thus conclude that the Hamiltonian describing the normal-metal is equivalent, in 
each ((fc, j'), (—fc, 4-)) subspace, to the Hamiltonian of a two-level system under the action 
of a circularly polarised held. Moreover, since TLs is time-independent and commutes 
with TLnit), it does not contribute to (A.22). We can thus write this equation as 

Ql3iv,T) = Tr 


0) ^ (r, 0) 


(A.33) 


where 


^ (f, 0) = T exp 


-- / duUniu) 
h Jo 


(A.34) 


and Pnifd) = with Zn{(3) = Tr(e“^^"^°^). Therefore, we have 

shown that there exists an equivalence, in what concerns the calculation of the work 
characteristic-function in the adiabatic limit, between the dynamics of the NIS junction 
and that of an assembly of independent two-level systems, subjected to a circularly 
polarised held. 

In Appendix B we will consider the calculation of the work characteristic-function 
of the latter system. 
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Appendix B. The work characteristic-function of a two-level system in 
contact with a thermal bath 


We now wish to compute (A.33) where TLnit) is given by (A.32). Since the different k 


modes are independent, we can restrict the calculation of the said function to that of a 
single mode. Furthermore, we will drop the k label, as no confusion can arise. 

The derivation of [Appendix A shows that one can reduce the dynamics of the 
system under the applied voltage protocol to the dynamics within the low lying energy 
states of the normal metal if one assumes that the transition time r' 3> h/A, see eq. 
(A.32). We will now further assume that <F(t) ~ hujt/2, where u = 2el4aax/h. Such 


an approximation, if it were exact, would actually imply that the voltage was being 
quenched. For it to be valid, one must have eVmax'rV^ ^ 1) i-6- t' ^ h/A„, as 
eVmax ~ An- Therefore, the time evolution of the system is only adiabatic with respect 
to the larger energy gap, being diabatic |102j with respect to the much smaller energy 
gap, i.e. one must have h/A -C r' -C h/A” ^ r. Therefore, for a clear separation 
between the two time-scales associated with the dynamics of the system to exist, one 
must have t'^(k) A, as stated above. 

We will change somewhat the notation with regard to the previous section, so as 
to keep the result obtained as general as possible, rather than identifying it solely with 
the dynamics of the NIS junction. In the new notation, the initial Hamiltonian is given 
by 


HniO) = ^ S' + ^ S' 


(B.l) 


where in the NIS junction context h = 2^” and F = 2A"', and where, without loss of 
generality. A"’ can be chosen to be a real-number, as one can always £x the arbitrary 
superconducting phase ip to be zero. 

The quantities h and F can be viewed as the components of a constant (pseudo) 
magnetic field applied to the two-level system. Since the system is initially in 
equilibrium with a thermal bath, the system’s partition function is given by Zn{fd) = 
2 cosh (/?ynT7?/2). 

The time-dependent Hamiltonian TLnit) can be written as 


^ h T 

Tinit) = 2^^ ~2 ^ cos{ut) + sin(a;t) ) 


(B,2) 


The dynamics of the two-level system is such that the constant applied held is 
substituted by a circularly polarised held in [0,r], i.e. by a Rabi dynamics during 
that interval. After the application of the protocol, the circularly polarised field is 
again replaced by a time-independent field and the system is described by the constant 
Hamiltonian In the process, the applied held has been rotated by an angle 

6 = ojT around the 2 ; axis. 

It is well known that for the Rabi dynamics, one can write the time-evolution 
operator ^(r, 0) in (A.33) as ^(r, 0) = with IZr = A representing a 
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pure rotation around the z axis and 'Wr = e where TL' is a time-independent 

pseudo-Hamiltonian, which depends on cn, and is given by 




(B.3) 


Substituting this expression for ^ (r, 0) in (A.33) and noting that Tl\ TZr = 

^iTzt-HniT)Rru _ since = 'Hn(O), one obtains for the characteristic- 

function the result 




Znifd) 


Tr 




(B.4) 


Performing the trace over the complete set of states that diagonalises 'Hn(O), we obtain 
for ^^(u,r), the expression 

1 


G/siv^r) = 


Znifd) tf, 


E I (n 


jn'r/hl / 


a',n) 


X g*(f^'-f^)Wr2+h2/2g-/3tTv/r2+h2/2 


(B.5) 


where the unit vector n refers to the direction of the applied held at f = 0, its 
components being given by fix = Moreover, writing 

n' = |7r2 + (h-hcu)2(n' ■ a), where n' = . ^ and n' = . , we 

can expand the exponential as 


eiWrin ^ 


2h 


T^Jv^ + {h — huy 

W 


+ i[\\ ■ fi) sin 

+ [h — hu) 


n • cr 


i sm 


2h 


[fi X (n' X n)] ■ CT . 


(B.6) 


The hrst two terms of (B.6) are diagonal in the | cr, n) basis, whereas the last term is 
only non-zero when evaluated between two states of the | cr, n) basis with opposite spin. 
We thus obtain that the square moduli of the matrix elements that appear in eq. (B.5) 
are given by 




2 _ 

2 / 


COS 


ry -|- (h — huY 

“ 2h ”” 


fn' ■ n)^ sin^ 


ry'r2 + (h - huf 
2h 


I n X (fi' X fi) 1^ 

r-y/r^ + {h — huY 


X sm 


2h 


drr-o 


(B.7) 
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Taking into account that | n x (fi' x fi) p= 1 — (n' ■ fi)^ and using the expressions for 


the components of n and ii' given above in (B.7), we dually obtain the expression 

-jfrpih 


Q0{v,t) = 


cos 


hooy 


2h 


X 


(r^ + - hhuf 

(r2 + h2)(r2 + 

(r2 + h2)(r2 + (h-hn;)2) 

cosh {(3/2 + iv) A/r2 + ^2 


■ sin 


sin 


+ {h — huy 
^r2 + (h-hn;)2' 


2h 


cosh (/?Vr2 + h2/2) 


(B.8) 


It is trivial to check that Gpiv, r), as given by (B.8), fuldls both the equality r) = 1 
and Gp{i(3^T) = 1, which provides a check on the correctness of the result, since these 
equalities are built into the dednition of r) by construction, as pointed out above. 

Performing the inverse Fourier transform on Qp{u,T), we obtain for P^(IF, r) the 
result 


F^(IP,r) = 


cos 


yr2 + (h-hn;)2' 
2h 


(F^ + - hhuf 


(F2 + h2)(F2 + (h- 
F2(hn;)2 


hen) 2 


■ sm 


(F2 + h2)(F2 + (h-hcn)2) 

P 0 W /2 


sm 


ry'F2 + {h- huf 
2h 

v'r" + (A 


5(hF) 


r 




2h 


2 cosh (/3v/r2 + h2/2) 

5 (W - VF2 + h2) + 5 (it + VF2 + h2) 


(B.9) 


The form of (B.9) is easily interpreted from the two-level structure of the problem, as 


either no work is performed on the system if the time-dependent held does not induce 
transitions between the levels (drst term), or otherwise the time-dependent held induces 
a transition between the ground state and the excited state, involving an amount of work 
performed on the system equal to IT = V^F^TpTp, or the inverse transition is induced 
involving a negative amount of work IT = — \/F2 -|- being performed on the system 


(second term). Note that if one sets IT —)■ —IT, u —>■ —cu and h —)■ —h, F —>■ —F in (B.9), 


one can directly check that P/s{—W,—t) = e~^^ PjsfW^r), i.e. the work distribution 
satisdes the Crooks-Tasaki relation. 

The average work performed on the system during the application of the protocol 


can be computed either by diderentiating (B.8) with respect to u, or directly from (B.9), 
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T^ihui) 


sin 


+ {h-hu:Y' 


2h 


VT^TW{r^ + {h - huf) 

X tanh (/3 VpT7?/2) . (B.IO) 

Since the change in the free energy of the system is zero dnring the transformation, this 
quantity is equal to the energy dissipated by the system into the thermal bath during 
the equilibration process occurring after t > t. One can easily check that, according 
to ( B.10[ ), {W)p is always larger or equal to zero in agreement with the second law of 
thermodynamics. It has a maximum at resonance, i.e. if ca = (F^ + h‘^)/{hh) and if 
car I r I /\/r2~+7? = (2n + l)7r, with n being an integer. Note that for a single two-level 
system, one can, for each value of r, choose oj such that (hF )/3 = 0, i.e. one takes u such 
that r^/r^ [h — huY/% = 2nn. The existence of such a minimum can be understood 

from that fact that for such choice of ca, = 1, and thus no transitions between 

the levels are induced by the unitary transformation, see also (B.5). This behaviour 


regarding the dissipated work in a two-level system mirrors the corresponding behaviour 
of the Rabi formula for the transition probabilities of such a system under the influence 
of a circularly polarised held. 

The mean-square deviation of the work performed during the application of the 


protocol can be computed either by differentiating (B.8) twice with respect to n, or 


directly from (B.9), and is given by 


{mY)0 = 


(r2 + (h - huf) 


sm 


Yr2 + {h- huf 
2h 


X 


1 - 


T^{nuf 


sm 


rwr2 {h — huY 


2h 


(B.ll) 


(r2 + h2)(r2 + (h-hc^)2) 

X tanh (/3\/r2 + h2/2)' . 

In Sect. 1^ we use (B.8) to compute the work dissipated by the NIS junction 
due to the application of the voltage protocol. Since it is also of interest, albeit not 


for the solution of the NIS problem, we leave to Appendix C the calculation of the 
work characteristic-function for an isolated two-level system due to the application of a 
circularly polarised held to such a system. 


Appendix C. The generating function of an isolated two-level system 


One can use the results obtained in Appendix B to compute the characteristic function 
or the work distribution function for an isolated system, described by a micro-canonical 
ensemble, which undergoes a transformation that is analogous to the one described in 
that section, i.e. the system is initially isolated and is coupled to the circularly polarised 
held at t = 0, being decoupled from it at f = r. Such characteristic function is given 
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by an inverse Laplace transform |19] that involves G/ 3 {v,t) and the partition fnnction 
Z{P) = 2 cosh (/?Vr2 + h^/2) 

t) u:o{E) = ^ ^ > (C-1) 

where E is the energy of the system and ojq{E) = 6 (^E — a/L^ + h‘^/2^ + 

6[e + v/PT 7?/2) is the density of states of the system at f = 0. The contour C 
is the inverse Laplace transform contour from c — ioo and c + ioo, with c chosen such 
that all the singularities of the integrand are located to the left of c. In our case, c = 0. 

Note that these results cannot be generalised to an assembly of two-level 
systems, since the work characteristic-function is in this case the product of 

characteristic-functions for the individual systems. 

The integral can be readily performed and after factoring the term u:q{E) out, we 
obtain for t) the result 


A Jr^T^ + ih-huy 

fo(^..r)= - - - 

(T^ + - hhuf 

(r2 + h2)(r2 + (h-hn;)2) 


+ (C.2) 

ryr2 + {h- huf 


sm 


sm 


2h 

V'y/r 2 + [h — huy 


+ 


—2ivE 


2h 


(r2 + h2)(r2 + (h-hn;)2) 

From this expression, one can obtain, as above, by inverse Fourier transformation of 
Ge{.v,t), the work function distribution for an isolated two-level system. This is given 
by 


Pe{W,t) = 


cos 


+ 


sm 


ry F2 + (h — huY 

2h 

(F2 + - hhuf 

(F2 + h2)(F2 + (h-hn;)2) 

T^ihuf 

^ (F2 + h2)(F2 + (h-hn;)2) 

+ {h-hu:Y' 


+ {h- huf 
2h 


5{W) 


X sm 


2h 


5{W + 2E) 


(C.3) 


This result can be easily interpreted if one again notes that either no transitions occur 
between the two levels and in this case the work performed is zero (hrst term), or 
otherwise the work performed is —2E where E is the energy of the initial level (second 
term). Also, note that Pe+w{.—W^ —t) = Pe(IF, r), which is the version of the Crooks- 
Tasaki relation appropriate for isolated systems in which the expressions for the micro- 
canonical density of states corresponding to the initial and to the hnal Hamiltonian are 
identical, since the spectrum of the system does not change under application of the 
work protocol [H]. 
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We can also compute the average work dissipated in the transformation, using the 
expression for t) or that for Pe^W, t). This is given by 


{W)e = 


2ET^{huj) 


(r2 + h2)(r2 + (h-hn;)2) 


sm 


Vy^r2 + (h — huy 


2h 


(C.4) 


Again, one can minimise or maximise this quantity by an appropriate choice of the value 
of uj. This result may be relevant in the context of quantum phase-shift gates [103] where 
the transformation can be implemented at hnite frequency without generation of heat. 

Finally, we can also compute the mean-square deviation of the work performed 
during the application of the protocol, using the expression for Qe{,v,t) or that for 
PEiW, t). This quantity is given by 


{mr)E = 


AE^T^{huj) 


(r2 + h2)(r2 + (h-hcn)2; 


■ sm 


'+ {h - huy 

"W 


(C.5) 


1 - 


(r2 + h2)(r2 + (h-hn;)2) 


sm 


y^r2 + {h- huf 
2h 


Appendix D. tunnelling matrix element extraction from experiment 


The power dissipated per atom in a NIS junction in the steady state, when submitted 
to a constant voltage, is given by [HU ESI ESI EH EH] 

r = l‘”(k)P(/(S - eV) - /(®) - eV) 

at 

+ \v^yS{ftf + Et,-eV)], 

where /(e) = -I- 1) is the Fermi function and where 


(D.l) 


|MkP =1(1 + 4 


= 4(1-4 


(D.2) 


are the squares of the occupation factors in the superconductor. Using the density of 
states per atom (and per spin) in the normal metal p{f) introduced above, we can write 
this equation as 


V = 


AttcV r+°° 


df 






Ve + A2(0 

where the function 7 (,^) is given by 

7(0 = (v+++0 + eV/2) -eV- + A^K)) 

+ (y+T+O - er/2) iK -eV+ v+Fa+)) . 


(D,3) 


(D,4) 
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Assuming that the gap in the superconductor is constant within the Debye window 
of frequencies, one can perform the above integral, where only the second delta function 
in the dehnition of 7 (^) gives a contribution. We obtain 

27r |f(Co)P p(^o) sinh(e|l/|/2fcBT) 


V = 


(D.5) 


he|l/|[cosh(e|l/|/2fcBT) + cosh(AV2e|l/|fcBT)] ’ 
where .^o = &V/2 — 12^ j2eV. This implies that in hrst-order perturbation theory there 
is no power dissipated, and hence no current flowing for voltages well below the gap. 
Note, however, that eq. (15) corresponds to a second-order calculation, since A"" oc \t\^. 

Note that in the model described by eq. (|^, the tunnelling process conserves 
momentum in the plane of the dims, and hence the junction is not ohmic when both 
these dims are in the normal state. However, one can measure the power dissipated 
when e|H| = A in the superconducting state of the aluminium dim. We obtain from 

(ra 


V = ^ |f(0) I V(0) A tanh(A/2A;sT). (D.6) 

Therefore, a measurement of V at e|H| = A, as a function of the temperature, allows 
the extraction of |t(0)p. Note that A(T) is itself a function of the temperature, see [92]. 

Note that an expression of the tunnelling matrix element in terms of parameters 
of the system such as the thickness of the insulating layer (see comment on section 
always requires a model of the potential barrier, see e.g. una. The above method avoids 
the need to resort to such modelling. 


Appendix E. Techniques for measurement of heat released in the NIS 
junction 


Several techniques can be used to measure the small heat release that was computed in 
section The most straightforward technique would be the use of standard calorimetry 
[num. In such a case, if the conversion of the work performed on the dim, into heat, 
occurs in a short time frame after the end of the work protocol, there is an increase AT 
in the temperature of the normal metal dim given by 


AT = 




(E.l) 


ncy 

where n is the number of moles of material contained in the probe, and Cy the material’s 
molar heat capacity. If we use (W)p ~ 0.7nJ/g, as computed in section]^ we obtain a 
variation of temperature AT = 1.28 x 10“^ K, likely to be too small to be measured. 

A possible alternative technique would be the measurement of the released heat 
at constant temperature by coupling the system to an external probe that absorbs a 
quantity of heat, AQ = (W)p, isothermally. This transformation would be analogous to 
the isothermal expansion phase of a Stirling engine |1U5] . If an indnitesimal amount of 
heat 5Q is released in the NIS junction during the equilibration process, after application 
of the protocol (as stated in section the equilibration time is of the order of the 
electron-phonon relaxation time Tpe), we have from the second law of thermodynamics. 
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that 5Q/T = dS-pr + dS^iQ, where dSpr is the variation of the entropy of the probe 
(which would be the elastic degrees of freedom of the normal metal him) and diSHB is the 
variation of entropy of the heat bath (surrounding helium bath). Since the temperature, 
volume and number of atoms of the bath do not change, dS'jjp = 0. Thus, we have 


5Q/T = dSpr = 


dS 

Wr 


dy = 


dP 

Wr 


dr = 


r 


ar 


dT 

\p 

ar 

ap 

T 


dr = ^dr, (E.2) 

kj' 


where we have applied the Maxwell relation on going from the second to the third 
equality and the implicit function theorem on going from the third to the fourth equality 
in the equation above, and where kt is the isothermal compressibility of the probe and 
a its thermal expansion coefficient, as follows from equilibrium thermodynamics. The 
reasoning above is analogous to the one made when discussing the measurement of the 
magnetocaloric effect from magnetisation data in isothermal conditions, see Ref. [106] . 
Note that the stresses in the him and the pressure of the fluid that makes up the heat 
bath do not need to be equal and moreover, that the mechanical element that controls 
the expansion of the him would only move in time scales much larger than Tpe, necessary 
for the equilibration of the elastic degrees of freedom of the him. 

Assuming that the thermodynamic quantities that enter eq. (E.2) are independent 
of the volume for the small variations measured, we obtain, integrating this equation 

Kt 


Ar = ^AQ. 
al 


(E.3) 


Note that in this case the variation of volume is inversely proportional to the thermal 
expansion parameter, which implies that the smallness of this quantity in aluminium at 
low temperatures will actually amplify the ehect to be measured. 

Therefore, the (average) relative variation of the volume of the probe is given by 

(E.4, 

r 'ynCyT 

where 7 = is the Gruneisen parameter of the material that constitutes the probe 
and Vm is its molar volume. For aluminium, 7 1.7 ma and c„ ~ 1.5mJmol ^ 


The measurement of the relative variation of the probe’s volume in isothermal 
conditions thus gives direct experimental access to the work produced during the 
prescribed protocol. A histogram of the distribution of data obtained for multiple 
realizations of the protocol will allow the reconstruction of the full work distribution. 

If we use {W)i 3 ~ 0.7nJ/g, we obtain a relative deviation ^ 7.5 x 10“®, 

on average, which should be measurable using, e.g. capacitive methods or SQUID 
dilatometers, as such a relative deviation is in the limit of precision of capacitive methods 
[niiiH], provided that other sources of heat dissipation can be minimized. 
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